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I.  INTRODUCTION 


Ozone  has  intense  absorption  in  the  middle  ultraviolet  region  and 
serves  to  shield  the  earth  from  the  sun's  radiation  in  this  region.  A 
knowledge  of  the  spacial  and  temporal  distribution  of  ozone  is  important 
in  the  study  of  the  chemistry  of  the  atmosphere  and  in  the  evaluation  of 
the  impact  of  man-made  catalytic  trace  constituents  upon  the  ozone 
shield. 

The  attenuation  of  solar  energy  by  ozone  may  be  used  to  determine 
the  height  distribution  of  ozone  in  the  upper  atmosphere.  The  optical 
measurements  needed  to  determine  the  attenuation  are  made  by  a sequential- 
sampling  ultraviolet  filter  photometer  which  is  carried  aloft  by  a 
rocket  and  which  descends  on  a parachute.  The  design  of  the  ozonesonde 
has  been  described; 1 the  system  being  used  currently  is  a modification* 
of  that  described  in  the  reference.  The  ozone  concentration  measure- 
ments are  made  at  altitudes  of  20  to  70  km. 

This  report  describes  the  measurement  of  the  transmittance  of 
filters  used  in  the  rocketborne  ozonesonde  and  the  treatment  of  the  data 
by  computer  to  calculate  effective  ozone  absorption  coefficients.  The 
application  and  operation  of  the  programs  are  described.  The  report  is 
intended  for  those  interested  in  the  measurement  procedures  and  calcula- 
tions used  in  obtaining  effective  ozone  coefficients.  Program  listings 
are  not  given  but  may  be  obtained  upon  request. 


II.  EFFECTIVE  OZONE  ABSORPTION  COEFFICIENT 


O 

The  maximum  ozone  absorption  occurs  a£  2553  A,  and  measurable 
absorption  is  found  to  approximately  3300  A.  If  the  attenuation  of 
direct  sunlight  by  the  earth's  atmosphere  is  assumed  to  be  due  solely 
to  absorption  by  ozone  and  to  Rayleigh  scattering,  then  the  intensity 
of  direct  sunlight  is  given  by  the  expression  (Beer's  law) 


The  modified  system  will  be  described  in  a NASA  report  by  Krueger 

et  al. 

1 Naval  Weapons  Center.  Rocket  Ozonesonde  (ROCOZ)  — Design  ind 
Deoelojment , by  A.  J.  Krueger  and  W.  R.  McBride.  China  Lake,  Calif., 
NWC,  July  1968.  40  pp.  (NWC  TP  4512,  publication  UNCLASSIFIED.) 
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K*. h) 


t co  -u(A)*x(h)  -0(A)m(h) 
I ( A ) e e 


(1) 


where  I(A,h)  is  the  intensity  at  wavelength  A and  height  h (the  other 
quantities  are  defined  below).  Since  the  filters  have  finite  optical 
bandwidths,  the  monochromatic  absorption  coefficients  a(A)  are  replaced 
by  an  effective  absorption  coefficient  which  is  defined  as  the  weighted 
mean  value  of  the  ozone  absorption  coefficients. 


a . (x,m)  In 

3 x 


Si  (A)Q(A)S(A)F.(A)e"a(X)'VtS(X)md> 
2 J 

/I  (A)Q(A)S(A)F. (A)e_0(X)mdA 
o 3 


(2) 


where 


ctj  (x,m) 

j 

x 

m 

A 

I0(A) 
Fj  (A) 
Q(A) 
S(A) 
a(A) 
6(A) 


= effective  ozone  absorption  coefficient  for  the  j*”*1  filter 
= filter  number 

= slant  path  ozone  amount  from  height  h above  the  earth  to 
the  top  of  the  atmosphere 

= slant  path  air  mass  above  height  h 

= wavelength 

= extraterrestrial  solar  irradiance 
= transmittance  function  for  the  j ^ filter 
= photodetector  response 

= transmittance  function  for  the  diffuser  plate 
= ozone  absorption  coefficient 
= Rayleigh  scattering  coefficient. 


This  expression  is  evaluated  numerically  at  1-A  intervals. 
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The  values  employed  for  the  extraterrestrial  solar  irradiance,  the 
ozone  absorption  coefficients,  the  Rayleigh  scattering  coefficients,  and 
the  photodetector  response  were  supplied  by  A.  J.  Krueger*  of  NASA, 
Goddard  Space  Flight  Center,  Greenbelt,  Maryland.  The  transmittance 
function  of  the  diffuser  S(a)  and  the  transmittance  function  of  the 
filters  F(A)  come  from  laboratory  measurements  of  all  the  filters  and  a 
few  of  the  diffusers  as  described  in  the  next  section. 


III.  TRANSMITTANCE  MEASUREMENTS 


The  transmittance  of  each  filter  and  some  of  the  diffusers  is 
measured  with  a dual  range  (0  to  100%,  0 to  10%)  Cary  Model  14R  spectro- 
photometer equipped  with  two  encoders  which  convert  shaft  positions  to 
digital  form.  Both  analog  and  digital  records  of  the  data  are  made. 

The  pen  position  encoder  (transmittance)  is  coupled  to  the  pen 
drive  mechanism;  the  encoder  goes  from  -50  to  +1050  counts  for  360°  of 
shaft  rotation.  A span  of  1000  counts  corresponds  to  a range  of  0 to 
100%  transmittance  (T)  or  0 to  10%T  depending  upon  the  position  of  the 
range  selection  switch  of  the  spectrophotometer.  Thus,  a change  of  one 
count  corresponds  to  a change  of  0.1%T  on  the  0 to  100%T  range  and  to 
0.01%T  on  the  0 to  10%T  range.  The  zero  of  the  pen  (and  hence,  of  the 
analog  record)  seldom  corresponds  to  the  zero  of  the  encoder.  For  exam- 
ple, the  encoder  may  give  28  counts  when  the  pen  is  at  its  zero  position. 
Then,  on  the  0 to  100%T  range,  an  encoder  reading  of  1028  counts  corre- 
sponds to  100%T.  This  offset  between  the  zeroes  of  the  pen  and  the 
encoder  must  be  corrected  for;  the  procedure  and  program  for  doing  this 
are  described  later  in  this  report.  The  wavelength  mechanism  of  the 


The  data  supplied  were  arrived  at  in  the  following  manner.  The 
solar  irradiance  values  were  obtained  by  averaging  a high  resolution 
spectrum  (Furukawa  et  al.)  over  1-A  intervals.  The  photodetector 
response  was  taken  as  a constant  and  set  equal  to  one.  The  ozone  absorp- 
tion coefficients  were  taken  from  Vigroux  and  Inn  and  Tanaka. '*  The 
Rayleigh  scattering  coefficients  were  calculated. 1 

National  Center  for  Atmospheric  Research.  Composite , High  Reso- 
lution Solar  Spectrum  from  2080  to  3600  A,  by  P.  M.  Furukawa,  P.  L. 
Haagenson,  and  M.  J.  Scharberg.  Boulder,  Colo.,  NCAR,  February  1967. 
(NCAR  Technical  Note  No.  26,  publication  UNCLASSIFIED.) 

E.  Vigroux.  "Contribution  a H'etude  exper imentale  de  ''absorp- 
tion de  '’.'ozone,"  Ann.  Phys. , Vol.  8 (1953),  pp.  709-62. 

Edward  C.  Y.  Inn  and  Yoshio  Tanaka.  "Ozone  Absorption  Coeffi- 
cients in  the  Visible  and  Ultraviolet  Regions,"  in  Ozone  Chemistry  and 
Technology.  Washington,  D.C.,  American  Chemical  Society,  1959.  Advatices 
in  'hemietry  Serine  >C:],  pp.  263-68. 
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spectrophotometer  is  connected  directly  to  an  encoder  which  converts 
shaft  position  to  digital  form.  One  count  on  the  encoder  corresponds 
to  1 A.  The  output  of  the  encoders  goes  to  a Datex  SDS-1  spectropho- 
tometer data  recording  system  having  a punched  paper  tape  output.  The 
control  module  of  the  SDS-1  system  provides  for  the  command  to  read 
out  data  and  for  the  recording  on  tape  of  identification  information  in 
the  form  of  a 10-digit  number,  each  digit  of  which  is  manually  set  on 
one  of  10  thumbwheel  parameter  switches.  The  system  uses  the  IBM  odd 
parity  code  given  in  Table  1 of  Appendix  A. 

The  10-digit  identification  number  is  recorded  twice  at  the  begin- 
ning of  every  spectrum  recorded  on  tape.  The  first  six  digits  give  the 
date  the  spectrum  was  recorded  (month,  day,  year),  the  next  two  digits 
designate  the  filter  assembly,  and  the  last  two  digits  give  the  spectrum 
number.  Thus  the  identification  number  1024751206  identifies  the  record 
which  follows  as  being  the  6th  spectrum  run  on  24  October  1975  and  shows 
that  it  is  for  filter  assembly  number  12.  Next  the  transmittance  (in 
digital  counts)  and  wavelength  data  are  recorded  every  angstrom.  Tables 
showing  the  identification  number  record  format  and  the  transmittance- 
wavelength  data  record  format  are  given  in  Appendix  A. 

A.  TRANSMITTANCE  OF  FILTERS* 

Narrow-band  optical  filters  having  bandpasses  centered  at  approxi- 
mately 2630,  2880,  3010,  and  3200  A and  having  half-widths  of  about 
32  to  41  A except  for  the  2630  A narrow-band  filter  which  has  a half- 
width  of  about  103  to  109  A,  are  used  in  conjunction  with  a broad-band 
filter  to  isolate  narrow  regions  of  the  spectrum.  The  broad-band  filter 
is  made  of  single  crystalline  nickel  sulfate  hexahydrate  sandwiched 
between  plates  of  Corning  9863  glass  and  is  used  to  isolate  the  ultra- 
violet region  and  block  the  visible  region.  The  narrow-band  filters  are 
mounted  in  a rotatable  filter  wheel.  The  filter  wheel  and  broad-band 
filter  are  mounted  in  the  sample  compartment  of  the  spectrophotometer  in 
such  a manner  that  the  monochromatic  radiation  passes  through  the  narrow- 
band  filter  (mirror  side  towards  incoming  radiation),  a 3/8-inch  stop, 
and  then  the  broad-band  absorption  filter.  The  spectra  are  scanned  at 
2.5  A/sec  from  3425  to  2500  A.  A short  slit  height  is  used  to  reduce 
the  chance  of  inadvertent  beam  clipping  by  the  filter  wheel. 


* The  term  filter  without  a modifier  specifying  narrow  or  broad 
will  refer  to  a filter  composed  of  two  parts,  one  a narrow-band  filter 
and  one  a broad-band  filter.  Thus,  the  term  3200-A  filter  refers  to 
the  combination  of  a 3200  A narrow-band  filter  with  a broad-band 
filter. 
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The  long  wavelength  tail  of  the  2630-X  filter  has  a small  light 
leak  which  lies  at  a wavelength  at  which  the  atmosphere  is  far  more 
transparent  than  it  is  at  2630  A.  Hence,  the  leak  is  important  in  the 
calculation  of  ozone  concentrations.  In  order  to  examine  the  long 
wavelength  tail  of  the  2650-A  filter  on  an  extended  range  a neutral 
density  screen  of  about  1%T  is  used  to  attenuate  the  reference  beam  of 
the  spectrophotometer;  this  results  in  full  scale  on  the  0 to  10%  range 
of  the  instrument  being  approximately  0.1%T. 

The  transmittance  of  most  of  the  filters  was  measured  on  the  0 to 
10%T  range  except  for  a few  whose  transmittance  exceeded  10%.  Because 
the  100%T  line  of  the  instrument  could  not  be  made  completely  flat 
(i.e.,  independent  of  wavelength)  with  the  available  multipot  adjust- 
ments of  the  spectrophotometer,  the  100%  line  was  recorded  and  used  to 
correct  filter  spectra  for  small  deviations  of  the  100%  line  from  the 
value  100. 

B.  TRANSMITTANCE  OF  DIFFUSER 

The  integrating  sphere  of  the  original  design  of  the  ozonesonde  has 
been  replaced  with  a diffuser  plate.  Two  types  of  diffuser  assemblies 
have  been  used.  One  is  made  up  of  three  pieces  of  ground  fused  silica 
cemented  together.  The  other  is  a barium  sulfate  coated  fused  silica 
plate.  The  transmittance  of  the  diffusers  was  measured  on  the  0 to  10%T 
range  with  a neutral  density  screen  of  approximately  1%T  in  the  reference 
beam.  The  spectra  were  scanned  at  2.5  X/sec  from  3425  to  2500  X;  the 
full  slit  height  was  used.  Data  were  recorded  every  angstrom. 

IV.  CALCULATION  OF  FILTER  TRANSMITTANCE  FUNCTION  F(>) 

The  purpose  of  this  section  is  to  describe  the  methods  used  in  the 
computer  programs  to  transfer  data  from  a record  on  punched  paper  tape 
and  to  compute  numerical  quantities  characterizing  the  filter  transmit- 
tance, F . ( i ) . 

The  programs  were  coded  in  FORTRAN  V and  FLECS.*  The  plotting 
routines  all  use  the  Integrated  Graphics  System  and  the  FR-80  plotter  by 
means  of  the  Naval  Weapons  Center  (NWC)  subroutine  CURVEG  or  PLOTSG.' 


* FLECS  (Fortran  language  with  extended  central  structures).^1 

Terry  Beyer.  FLECS:  User's  Manual,  Computing  Center,  University 
of  Oregon. 

Naval  Weapons  Center.  User  Oriented  Computer  Graphics  With  Sub- 
POutine  PL0TSGt  CURVEG , and  SHOMAT,  by  Alan  L.  Craig  and  C.  Howard 
Shomnte.  China  Lake,  Calif.,  NWC,  March  1972.  90  pp.  (NWC  Technical 

Note  404-133,  publication  UNCLASSIFIED.) 
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The  programs  are  an  disc  storage  at  NWC  under  the  qualifier  "ROCOZ", 
with  file  names  TAPECARD,  EDIT,  and  OFFSET.  The  names  of  the  elements 
in  each  file  are  given  in  Appendix  B. 

The  process  of  obtaining  the  Fj  is  broken  into  several  stages  to 
permit  evaluation  and  convenient  editing  of  the  original  data.  This 
method  also  allows  modifications  for  observing  the  effects  of  variations. 
As  a result,  the  method  is  not  an  efficient  production  model,  nor  is  it 
intended  to  be. 

A complete  data  run  is  made  up  of  the  following  four  programs,  each 
of  which  will  be  discussed  separately:  the  TAPECARD  program  which  trans- 

fers spectrophotometric  data  from  punched  paper  tape  to  punched  cards; 
the  EDIT  program  which  checks  the  cards  for  bad  punches,  ordering  of 
wavelengths,  and  duplicate  or  missing  data;  the  OFFSET  program  which 
applies  an  offset  correction,  smooths  the  100%  line,  and  calculates 
percent  transmittance  (the  Fj's);  and  the  ALPHAEFF  program  which  com- 
putes the  effective  ozone  absorption  coefficient.  This  section  (IV) 
describes  the  first  three  computer  programs. 


A.  TAPE-TO-CARD  PROGRAM,  TAPECARD 

The  TAPECARD  program  converts  the  information  on  punched  paper  tape 
to  data  punched  on  cards.  The  tape  is  read  on  a PDP-11  computer  into  a 
named  storage  file  in  the  UNIVAC  1110  computer,  where  data  are  stored  in 
PDP  form.  Each  character  is  scanned  and  converted  from  paper  tape  code 
to  punched  card  code.  When  80  characters  have  been  scanned,  a line  is 
printed  and  a card  is  punched.  The  pattern  of  numbers  on  the  tape  is: 

hsxxlSxxxxx,  where  K is  a blank,  x is  a digit,  and  s is  either 
a digit  or  the  symbol  //. 

A full  line  of  printed  output  consists  of  eight  of  these  patterns.  If  a 
space  does  not  occur  where  it  should,  the  scan  for  that  80  characters 
stops  and  a new  scan  commences.  When  this  happens  the  punched  card  and 
printed  output  correspondingly  will  have  only  a partial  line,  followed 
by  a new  line.  In  this  way,  the  transmittance-wavelength  pairs  are 
aligned  making  these  errors  obvious  on  the  printed  page.  Subsequent 
card  corrections  are  also  easier  to  make. 

An  example  of  a card  deck  used  to  command  the  UNIVAC  1110  computer 
to  print  and  punch  the  data  received  from  the  PDP-11  computer  is  shown 
in  Figure  1.  The  only  changing  quantities  in  the  card  deck  are  the 
names  of  the  punched  paper  tapes.  The  tape  name  begins  in  column  1 and 
uses  no  more  than  six  characters.  The  tape  name  is  followed  by  a period 
and  TAP.  There  must  be  no  blanks  between  column  1 and  P of  TAP.  Example: 
For  a tape  named  MEH1,  the  tape  name  card  would  read 

MEH1.TAP  . 
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@RUN  026PCH  . . . 

@ASG , A 3130131*PDPPSEUD0001 . 
0ASG.AX  ROCOZ*TAPECARD . 

@XQT  ROCOZ*TAPECARD . PCH 
MEH1 . TAP 
MEH2 . TAP 


@FIN 

FIGURE  1.  Deck  Structure  for  the  TAPECARD  Program. 

B.  EDIT  PROGRAM 

The  EDIT  program  takes  raw  data  in  the  form  of  punched  cards  from 
the  TAPECARD  program,  edits  them  for  bad  punches  and  sequencing,  inter- 
prets coded  data,  and  tabulates  and  plots  the  data. 

Each  data  card  from  the  TAPECARD  program  contains  eight  transmittance- 
wavelength  pairs  in  fields  of  10.  A typical  card  for  the  100%  line  may 
contain,  for  example,  after  a first  blank,  the  sequence 

998  03363  #01  03362  #01  03361  999  03360  #03  03359 

#02  03358  #05  03357  #02  03356  #04  03355  #05  03354 

which  corresponds  to  transmittance  value  (in  digital  counts)  998  for 
wavelength  3363,  and  transmittance  1001  for  wavelength  3362.  The  "#" 
sign  is  a code  on  the  punched  paper  tape  which  is  changed  by  the  EDIT 
program  to  the  two  digits,  10.  The  program  searches  for  the  beginning 
wavelength  3399  X among  the  first  five  cards  (stopping  if  it  does  not 
find  the  wavelength).  Then  it  checks  whether  the  wavelengths  are  in 
descending  order  beginning  with  3399  X and  ending  at  or  before  2500  A. 

The  wavelength  and  its  corresponding  transmittance  value  are  stored  in 
arrays  IX  and  Y.  The  last  card  must  contain  a wavelength  _ 2500  A. 

While  the  program  checks  that  the  wavelengths  are  in  descending  order, 
it  also  checks  that  transmittance  values  (e.g.,  #13)  contain  the  proper 
characters  for  the  position,  and  whether  any  wavelengths  are  duplicated 
or  missing.  In  case  of  an  error,  an  appropriate  message  with  identify- 
ing wavelength  is  printed,  and  in  the  case  of  an  improperly  formatted 
transmittance  value,  a value  almost  equal  to  the  lower  limit  of  the 
plotting  scale  is  substituted,  so  that  the  error  shows  up  dramatically 
on  the  plot  (see  Figure  2).  If  twenty  errors  are  found,  the  program 
stops.  After  wavelength  2500  X has  been  reached,  the  wavelength- 
transmittance  value  pairs  are  printed  and  plotted. 
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FIGURE  2.  Plot  of  Raw  Data  by  the  EDIT  Program  for  the  J = 1 
Filter  of  the  LD-14  Filter  Assembly.  The  ordinate  scale  is 
approximate  because  the  zero  offset  and  the  100%  line  devia- 
tions have  not  been  corrected  for.  The  three  points  with 
zero  value  represent  bad  data. 

The  EDIT  program  treats  many  sets  of  data  identically.  Each  set 
requires  an  identifying  number  and  designation  of  the  spectrophotometer 
range  used  during  the  recording  of  data  in  order  for  the  program  to 
position  the  decimal  point  correctly  in  the  transmittance  data.  The 
identifying  number  of  the  spectrum  is  punched  out  by  the  TAPECARD  pro- 
gram at  the  beginning  of  the  transmittance  data. 

An  example  of  the  card  deck  used  to  run  the  EDIT  program  is  given 
in  Figure  3.  There  is  one  data  deck  for  each  spectrum;  as  many  data 
decks  as  desired  may  be  used. 
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The  cards  are  discussed  in  order  of  their  appearance. 
Card  1 ID  Card  Format  (6A6) 


Cols . 
2-11 


12-21 


Card  2 


Content 

ID,  identification  number  (e.g.,  0522751335;  this 
shows  that  the  data  were  taken  5/22/75  for  filter 
assembly  13,  spectrum  no.  35). 

ID,  identification  number  repeated  (not  required 
but  usually  present). 


IRANGE  Card 


Format  ( ) 


IRANGE,  LAMDEL,  NDEG  (e.g.,  3,  1,  0)  appearing 

anywhere  on  the  card  and  separated  by  two  commas. 

IRANGE  controls  the  placement  of  the  decimal  in  the  transmittance 
values  and  the  limits  of  the  ordinate  of  the  plot.  The  possible  values 
of  IRANGE  are  1,  2,  3,  or  4.  Table  1 shows  the  meaning  of  the  parameter. 


TABLE  1.  Values  of  IRANGE  for  the  EDIT  Program. 


Cary  14R  Ordinate  range 

1 KAiNv.iL  f , 

range  of  plot 


1 

0 

to 

100% 

-10 

to 

+110% 

2 

0 

to 

10% 

-l 

to 

+11% 

3 

0 

to 

100% 

80 

to 

120% 

4 

0 

to 

100% 

-5 

to 

+5% 

For  use  with 


spectrum  run  on  0 to  100%  range 
spectrum  run  on  0 to  10%  range 
100%  line  data 
0%  line  run  0 to'  100%  range 


LAMDEL  is  the  interval  in  angstroms  between  recorded  wavelengths.  This 
is  normally  1. 

NDEG  is  a plotting  parameter  either  left  blank,  set  to  zero,  or  set 
to  n,  in  which  case  a curve  of  the  nth  degree  is  fitted  to  the  data.  In 
this  case  the  parameter  M in  the  program  must  be  set  to  8 x NP  + 1000 
where  NP  is  the  number  of  points. 

Card  3 Data  Card  Format  (80R1) 

A data  card  contains  up  to  eight  transmittance- 
wavelength  pairs  punched  by  the  TAPECARD  program. 
Each  pair  occupies  ten  columns  in  the  pattern. 

tfsxxlixxxxx 

where  is  a blank, 

s is  either  a digit  or  the  symbol  It, 
and  x is  a digit. 
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The  output  of  the  EDIT  program  consists  of  a tabulation  of  a run- 
ning index  from  1 to  900  along  with  corresponding  wavelength  and  trans- 
mittance values.  Statements  giving  wavelengths  duplicated  or  missing 
are  printed  with  the  tabulation  as  are  statements  indicating  the  exist- 
ence of  nondigital  values  in  columns  which  should  contain  digital  values 
only.  A plot  of  transmittance  values  versus  wavelength  is  made  for  each 
spectrum  (Figure  2).  Erratic  values  show  plainly  on  the  graphical 
output . 

The  data  cards  are  corrected  as  indicated  by  the  EDIT  program  and 
are  used  to  calculate  the  transmittances  of  the  filters. 


C.  OFFSET  PROGRAM 

The  OFFSET  program  processes  seven  data  decks  to  obtain  the  percent 
transmittances  for  the  four  filters  making  up  a given  filter  assembly. 
The  first  two  programs,  TAPECARD  and  EDIT,  treat  all  the  data  the  same, 
merely  resulting  in  uncorrected  data  on  punched  paper  tape  being  changed 
into  corrected  data  on  cards.  In  the  OFFSET  program,  however,  the  data 
decks  all  belong  to  the  same  filter  assembly  and  must  be  arranged  (bar- 
ring exceptional  cases)  in  the  following  sequence. 

1.  100%  line  deck,  a reference  for  the  subsequent  decks 

2.  Nominal  3200  X filter  deck 

3.  Nominal  3000  X filter  deck 

A.  Nominal  2875  A filter  deck 

5.  Nominal  2650  X filter  deck 

6.  Nominal  2650  X filter  against  screen  B deck 

7.  Screen  B deck. 

The  expression  "filter  A against  screen  B"  means  that  the  transmittance 
of  filter  A was  measured  with  the  reference  beam  attenuated  by  a neutral 
density  screen  called  B.  Since  screen  B has  a transmittance  of  approx- 
imately 1%,  the  upper  limit  of  the  transmittance  scale  becomes  about 
0. 1%T. 


The  zero  of  the  spectrophotometer  pen  seldom  corresponds  to  the 
zero  of  the  pen  position  encoder.  The  difference  between  the  zeroes  of 
the  pen  and  the  encoder  is  called  the  offset  and  the  term  gives  rise  to 
the  name  of  the  program  which  corrects  for  it. 

The  OFFSET  program  subtracts  the  offset  correction  (an  input  quan- 
tity) from  each  point  of  the  100%  line,  smooths  the  line,  performs  the 
calculations  necessary  to  produce  percent  transmittances  for  each  filter 
band,  prints,  punches  cards,  and  plots  the  result.  The  following  out- 
line describes  in  algorithmic  form  the  instructions  the  computer  receives 
to  perform  the  necessary  calculations. 
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1.  Read  the  offset  correction  ZE100  and  the  wavelength-transmittance 
data  pairs  for  the  100%  line.  Call  the  EDIT  subroutine,  which  edits  the 
data  (see  the  description  for  the  second  program,  the  EDIT  program). 
Subtract  the  correction  ZE100  from  each  transmittance  value,  Y(I),  and 
scale  by  dividing  by  D = 10,  a value  which  is  assigned  by  IRANGE.  Find 
the  wavelength  nearest  the  2500  A end  which  is  a multiple  of  5 and  make 

it  the  last  wavelength.  Then  count  the  number  of  wavelength-transmittance 
pairs,  calling  the  total  N100. 

2.  Call  the  AVRIJ  subroutine  to  smooth  the  100%  line  by  making  a 
sliding  average  of  21  points,  assigning  the  average  value  to  the  mid- 
point. The  processed  values  are  called  SM100(I),  1=1,  N100. 

(The  program  then  moves  on  to  obtain  transmittance  values  correspond- 
ing to  data  in  six  decks  identified  by  JFIL  = 1,  •••,  6.  The  processing 
outlined  in  Step  3 is  identical  for  the  first  three  decks  but  the  last 
three  decks  are  treated  differently,  as  outlined  in  Steps  5,6,  and  7.) 

3.  Read  the  data  for  the  filter  consisting  of  JFIL  (filter  number) 
and  NZER.  NZER  is  the  input  offset  correction  for  the  filter  and  is 
left  blank  on  input  for  filters  JFIL  = 1,  2,  3,  4,  since  it  is  computed 
in  the  program  (see  below).  For  JFIL  = 5,  6,  a nonzero  offset  may  be 
input,  which  is  used  instead  of  the  computed  value. 


a.  Find  the  wavelength,  LMAX,  corresponding  to  the  maximum 
transmittance,  YMAX.  If  there  are  several  LMAXs,  choose  the  one  nearest 

O 

2500  A.  The  corresponding  index  is  called  IMAX. 

b.  Obtain  the  offset  correction  ZER  near  the  red  end  of  the 

O 

spectrum,  i.e.,  near  3400  A,  for  J = 1 through  5.  If  the  offset  correc- 
tion is  input  (i.e.,  NZER  ^ 0),  do  not  compute  ZER.  If  the  offset 
correction  is  not  input  (NZER  = 0),  average  the  50  transmittance  data 
corresponding  to  50  wavelengths  beginning  at  the  red  (long  wavelength) 
end  of  the  spectrum,  and  call  the  integer  closest  to  the  average,  ZER. 

c.  Average  the  transmittance  data  corresponding  to  the  last 
50  wavelengths  at  the  blue  end  of  the  spectrum,  i.e.,  near  2500  A,  and 
call  the  integer  closest  to  the  average  value,  ZEB,  the  offset  correc- 
tion for  the  blue  (short  wavelength)  end  of  the  spectrum. 


d.  If  ZER  = ZEB,  the  transmittance  TAU^(I)  is  given  by 
- (FJ(O-ZER)IOO 


TAU1(I) 


SM100(T) 


1=1,  NMIN 


(3) 


where  FJ(I)  is  the  transmittance  value  for  the  corresponding  wavelength 
LAM(T),  SM100(1)  is  the  value  of  the  smoothed  100%  line  corresponding  to 
wavelength  LAM(I),  NMIN  is  the  smallest  of  the  numbers  900,  N100  (the 
number  of  points  in  the  100%  line),  and  NN  (the  number  of  points  in  the 
filter  data) . 
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e.  If  ZER  4 ZEB,  subtract  ZER  from  each  transmittance  value 
near  the  3399  A end  to  LMAX,  and  subtract  the  blue  offset  value  ZEB  from 
LMAX  to  the  blue  end.  That  is, 

TAU1(I) 

and 

TAU1(I) 

f.  Reorder  the  points  so  that  the  first  wavelength 
LAM(l)  = 2500. 

g.  Find  the  maximum  transmittance,  YMAX,  for  the  corrected 
%Ts.  Determine  the  wavelengths  CLAMB  and  CLAMR  at  which  the  transmit- 
tance is  one-half  the  maximum  value,  YHAF  = YMAX/2.  Calculate  the 
center  of  the  bandpass,  CLAMC. 

CLAMC  = (CLAMB  + CLAMR) /2  . (6) 

Calculate  the  width,  W,  of  the  band  at  half-height. 

W = CLAMR  - CLAMB  . (7) 

h.  Punch  cards  with  wavelength  and  percent  transmittance 
pairs  for  the  four  filters  or  for  (2650/screen  B)  x screen  B in  a 
format  acceptable  to  the  ALPHAEFF  program,  along  with  the  constants  to 
be  used  in  that  program. 

4.  For  JFIL  = 2 and  3,  treat  the  3000  A filter  and  the  2875  A 
filter  as  in  3. 

5.  For  JFIL  = 4,  the  2650  A filter,  use  3a,  b,  d,  f,  g and  h.  Set 
ZEB  equal  to  ZER. 

6.  For  JFIL  = 5,  the  2650  A filter  against  screen  B,  (written 
2650/s  and  identified  as  JFIL  = 5)  has  only  the  red  end  of  the  band,  and 
hence  only  ZER.  Use  3a,  b,  d,  f and  g to  obtain  TAU^(I),  which  is 
stored  for  use  in  7. 

7.  For  JFIL  = 6,  for  screen  B,  set  ZER  = ZER4  if  NZER  = 0;  other- 
wise set  ZER  = NZER. 

Ca leu  Late 

TAU2(1)  - TAlipnl^^  . (8) 

The  desired  transmittance  for  the  long  wavelength  tail  of  the  2650  A 
filter  on  the  extended  range  is 

TAU(I)  = (TAU1(I)TAU^(I))100,  '(I)  > 2750  . (9) 


(FJ(Q-ZER)IOO 

SM100(I) 


I = 1 to  I 


MAX 


(4) 


(FJ (I)-ZEB) 100 
SM100(I) 


I = 


L1V+  1 to  NMIN.  (5) 
MAX 
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(If  either  TAU^(I)  or  TAU2(I)  is  nonexistent  for  a particular  I,  the 
product  is  set  to  zero  for  that  I.) 

The  input  data  decks,  seven  in  number,  are  the  corrected  decks  from 
the  EDIT  program.  The  input  quantities  are  those  of  the  EDIT  program 
plus  one  card  in  front  of  each  data  deck  to  notify  the  OFFSET  program 
what  data  deck  to  expect  (JFIL)  and  to  provide  an  offset  value  where 
necessary.  ZE100  is  the  offset  correction  for  the  100%  line.  It  is  a 
required  input  and  is  an  integer.  NZER  is  an  input  offset  for  2650/s 
and  screen  B decks.  It  is  not  required  and,  if  not  input,  the  offset 
will  be  computed  in  the  OFFSET  program.  It  should  be  noted  that  there 
is  a comma  after  ZE100  but  not  after  NZER  (see  Figure  4). 

An  illustration  of  the  input  decks  for  the  OFFSET  program  is  given 
in  Figure  4. 

Figures  5 through  8 show  the  transmittance  of  the  filters  for  LD-14 
as  calculated  and  plotted  by  the  OFFSET  program.  Figure  9 shows  the 
long  wavelength  tail  of  the  2650  A filter  on  an  expanded  scale. 


V.  CALCULATION  OF  DIFFUSER  TRANSMITTANCE  FUNCTION  S(A) 


The  purpose  of  this  section  is  to  describe  the  computer  program 
(called  DIFFUSER)  used  to  calculate  the  diffuser  transmittance  function 
from  data  stored  on  punched  cards  which  have  been  edited  by  means  of  the 
EDIT  program  described  in  Section  IVB. 

The  DIFFUSER  program  is  described  here  because  it  is  a modified 
OFFSET  program  and  in  its  logic  the  DIFFUSER  program  follows  the  OFFSET 
algorithm.  However,  there  are  only  three  input  data  decks  instead  of 
seven.  They  are  (1)  the  100%  line  deck;  (2)  diffuser  against  screen  B 
deck  (JFIL  = 5);  and  (3)  screen  B deck  (JFIL  = 6). 

The  input  quantities  are  ZE100  and  NZER  as  in  the  OFFSET  program. 

In  distinction  to  the  OFFSET  program,  NZER  is  a required  input  to  the 
DIFFUSER  program.  An  illustration  of  the  input  decks  for  the  program  is 
given  in  Figure  10. 

The  output  of  the  program  is  a tabulation  of  the  input  and  output 
data,  a plot  of  diffuser  transmittance  versus  wavelength,  and  a deck  of 
punched  cards  containing  the  diffuser  transmittance  in  the  format 
required  by  the  ALPHAEFF  program.  It  should  be  noted  that  although  the 
calculated  values  are  referred  to  as  transmittance  the  term  is  not  used 
in  its  usual  sense  as  the  spectrophotometer  is  designed  to  measure  the 
transmittance  of  nonscattering  samples  and  the  diffusers  are  highly 
scattering.  Figure  11  shows  the  transmittance  of  a BaS04  coated  type 
diffuser  and  Figure  12  shows  the  transmittance  of  a ground,  fused  silica 
type  diffuser.  The  DIFFUSER  file  is  listed  in  Appendix  B. 
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@RUN  0260FS  • • • 

@ASG,AX  K0C0Z*0FFSET . 

@USE  R,R0C0Z*0FFSET 
@XQT  R .MAIN 

C 100,26, 

0613751X010613751401 

3,1,0 

#21  03400  #21  03399  #23  03398  #23  03397  #21  03396  #21  03395  #24  03394  #19  03393 


#21  02503  #24  02502  #24  02501  #24  02500  #23  02499  #24  02498  #24  02497  #23  02496 
l^EOF 

'l. 

06137514030613751403 

2,1,0 

028  03402  028  03401  028  03400  028  03399  028  03398  028  03397  028  03396  028  03395 


^@E0F 

r2. 


JFIL  - 2 
(3000  A)  < 


JFIL  - 3 

(2880  A) 


JFIL  - 4 
(2630  A) 


JFIL  - 5 
(2630  A/S) 


JFIL  - 6 
(Screen) 


@EOF 

'3, 


@E0F 

r 

4, 


(3  EOF 
► 

5, 


@EOF 

6,29 


V?E0F 

$FIN 


FIGURE  4.  Deck  Structure  for  the  OFFSET  Program. 
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FIGURE  6.  Transmittance  of  the  J = 1 Filter  of 
the  LD-14  Filter  Assembly. 
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FIGURE  7.  Transmittance  of  the  J = 2 Filter  of 
the  LD-14  Filter  Assembly. 
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FIGURE  10.  Deck  Structure  for  the  DIFFUSER  Program. 
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FIGURE  12.  Transmittance  of  a Quartz 
Diffuser  Assembly. 
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VI.  CALCULATION  OF  EFFECTIVE  ABSORPTION  COEFFICIENT  i(x,m) 


The  effective  absorption  coefficient  ij(x,m)  for  filter  j is 
obtained  by  numerical  evaluation  of  the  expression 

. / 1 (V)Q(X)S(')F. (\)e~ ‘(')Xe_;  ( )mdA 

1,0  1 

- — In J - — , 

X /lo(^)Q(A)S(X)F.(A)e_r(  )md\ 

where  the  quantities  are  as  defined  in  Section  II  of  this  report.  The 
ALPHAEFF  program  carries  out  the  numerical  evaluation  and  then  approxi- 
mates the  effective  absorption  coefficient  by  a polynomial  as  described 
below. 


The  rectangular  integration  is  performed  by  multiplying  the  quanti- 
ties in  the  numerator  for  fixed  values  of  x,  and  m by  L (usually  1 A) 
and  summing  over  the  A’s;  the  computation  of  the  denominator  follows  the 
same  pattern.  The  two  sums  are  divided,  the  logarithm  (base  e)  of  the 
dividend  is  taken  and  then  multiplied  by  - 1/x,  yielding  ij(x,m)  for  a 
given  filter  j and  fixed  values  of  x and  m.  Whenever  x is  very  small, 
l'Hopital’s  rule  is  invoked  to  obtain  a limiting  form  of  the  formula  for 
ij(x,m).  The  values  of  uj(x,m)  are  computed  for  every  combination  of 
the  specified  values  of  x and  m (see  Appendix  C for  a table  of  x and  m 
values  used),  and  a least-squares  multiple-linear-regression  which 
minimizes  the  sums  of  the  squares  of  the  residuals  is  made  over  all 
computed  values  of  aj (x,m)  to  give  the  coefficients  of  the  polynomial 


2 2 

P.(x,m)  = a + a,x  + b,m  + a„x  + b„ra 
J o 1 1 2 2 


(10) 


For  the  initial  values  of  m (i.e.,  mG) , ;ij(x,mQ)  is  plotted  as  a function 
of  x;  the  values  of  P^(x,mo)  are  also  plotted  on  the  same  graph. 

Since  the  assembly  of  the  cards  and  data  decks  used  in  the  ALPHAEFF 
program  is  more  complex  than  for  the  other  programs,  the  process  will  be 
described  in  greater  detail  than  was  done  for  the  other  programs. 
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All  of  the  input  parameters  for  the  ALPHAEFF  program  are  punched 
out  by  the  OFFSET  program.  The  input  quantities  I0(A),  Q(A),  S(A), 
a(A),  and  8(A)  are  stored  on  a disc  and  called  up  for  use  by  the  @ADD 
control  cards.  Figure  13  is  an  illustration  of  an  assembled  ALPHAEFF 
deck. 

The  parameter  cards  punched  by  the  OFFSET  program  will  be  discussed 
in  the  order  of  their  appearance  in  the  assembled  ALPHAEFF  deck. 


Card  1 A-Parameter  Card 
Cols. 


11-15 


16-20 


Format  (4F5.0,  5X,  II) 

Content 

ALAMO,  the  beginning  wavelength, 
e.g. , 2500  A. 

ALAMF,  the  ending  wavelength,  e.g., 
3399  X. 

DELLAM,  the  interval  between  wave- 
lengths, usually  1. 

AK,  a scaling  factor  for  the  numerator 
integrand,  usually  1. 

NFLTRS,  the  number  of  filters  being 
processed,  usually  4. 


Card  2 Title  Card 
Cols. 


Format  (4A5,  15) 

Content 

Title,  e.g.,  LOKI  DART  12. 
Diffuser  number,  e.g.,  4. 


Card  3 x-Parameter  Card 


1-5 

6-15 

16-25 

26-30 

31-35 


26-40 


40-45 


Format  (15,  2F10.1,  F5.0,  F5.3,  15, 
F5.0) 

Content 

* 

J,  filter  number. 

CLAMC,  center  of  the  bandpass. 

W,  width  of  band  at  half-height. 

XINT,  initial  x value,  usually  0. 

DELX,  the  interval  between  x values 
(.05,  .02,  .005,  .001  for  filters 
0 to  3,  respectively,  are  in 
current  use) . 

NXS,  the  number  of  x's  to  be  computed 
(12,  24,  20,  20  currently  in  use). 

XFACT0,  a scaling  factor  for  the  x's 
(currently  XFACTO  = 1). 


* J,  the  filter  number.  On  input  J = 1,  2,  3,  4;  on  output  J = 0, 
1,  2,  3,  chiefly  for  historical  reasons. 
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Card  4 Identification  Card  Format  (10A6) 


Cols. 

2-11 

12-21 

Card  5 m-Parameter  Card 


10 

11-20 


24,25 


26-35 


Cols. 


Content 

ID  number,  e.g.,  0522751335. 

ID  number,  repeated. 

Format  (2F10.0,  15,  F10.0)  • 

Content 

m0,  initial  value  for  m,  usually  zero. 

DELM,  the  interval  between  computed 
values  of  m.  (The  values  .2,  .2, 
.01,  .01  are  currently  used  for 
filters  0 to  3,  respectively.) 

NMS,  the  number  of  m's  to  be  computed 
(10,  10,  05,  05  currently  in  use). 

FACTRM,  a scaling  factor  for  the  m's 
(currently  FACTRM  = 1). 


Card  6 MLR-Parameter  Card  Format  (513,  2F4.0) 


10-12  . 
13-15 

16-19 

20-23 


Content 

ITEST,  the  test  case  number  (0). 

N,  the  number  of  observations  (0). 

NP1,  the  number  of  independent  and 
dependent  variables  (5) . 

IW,  weighting  factor  indicator  (0). 

LIMIT,  the  maximum  number  of  steps 
in  the  regression  analysis  (5) . 

EFIN,  the  F value  for  entering  a 
variable  from  the  regression  (2.5) 

EFOUT,  the  F value  for  removing  a 
variable  from  the  regression  (2.5) 


Four  decks  from  the  OFFSET  program,  one  for  each  filter  of  the  filter 
assembly,  form  the  core  of  the  ALPHAEFF  deck. 

Except  for  two  beginning  and  two  ending  cards,  the  data  decks 
punched  by  the  OFFSET  program  consist  of  transmittance-wavelength  pairs 
punched  five  to  a card,  beginning  with  a wavelength  near  2500  A which  is 
a multiple  of  5 and  continuing  in  ascending  order  no  further  than  3399  X. 
These  cards  are  easily  distinguished  from  the  original  data  cards  (which 
were  in  descending  order  beginning  with  3399  X and  punched  eight  pairs 
to  a card) . 


* Values  of  the  parameter  currently  being  used  are  entered  in 
parentheses. 
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FICURE  13.  Deck  Structure  for  the  ALPHAEFF  Program. 
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If  noise  in  the  tails  of  the  transmittance  bands  was  not  removed 
previously,  it  is  removed  at  this  time  by  replacing  the  cards  having 
nonzero  transmittance  values  with  blank  cards  in  wavelength  ranges  where 
it  is  obvious  that  the  transmittance  should  be  zero.  Since  the  ALPHAEFF 
program  does  not  read  the  wavelengths,  blank  cards  may  be  used.  This 
does  away  with  the  necessity  to  punch  cards  and  for  this  reason  it  is 
easier  to  remove  noise  at  this  point  rather  than  earlier. 

Before  the  data  cards  are  used  in  the  ALPHAEFF  program,  it  is 
necessary  to  replace  part  of  the  cards  in  2650  A filter  deck  with  cards 
from  the  transmittance  measurements  made  with  a neutral  density  screen 
in  the  reference  beam  of  the  spectrophotometer.  The  latter  set  will  be 
called  the  "expanded  tail"  set.  The  crossover  from  normal  tail  to 
"expanded  tail"  is  made  near  0.09%T. 

An  end-of-file  card,  @E0F,  signals  the  end  of  the  data  for  a 
particular  filter. 

The  output  of  the  ALPHAEFF  program  consists  of  tabulations  and 
plots.  The  following  tabulations  are  produced. 

1.  List  of  the  input  data. 

2.  List  of  the  values  of  the  numerator  integrand  for  m = 0 and 
several  values  of  x (except  for  the  3200  A filter) . 

3.  Tabulation  of  x,tn  and  the  corresponding  values  of  aj(x,m). 

4.  Listing  of  the  input  parameters  to  the  multiple-linear- 
regression  program. 

5.  Tabulation  of  the  approximating  polynomial  coefficients  for 
each  iteration  until  a satisfactory  F-test  is  obtained. 

6.  Tabulation  of  Pj(x,m)  together  with  the  corresponding  values  of 
cij  (x,m) . 


Figure  14  shows  the  tabulation  of  approximating  polynomial  coefficients 
for  the  J = 1 filter  (nominal  3000  A filter)  of  LD-14. 

The  plots  are  of  two  kinds. 

1.  Plots  of  the  numerator  integrand  versus  wavelength  for  optical 
depths  approximately  equal  to  0,  1,  2,  3,  and  4 (except  for  the  3200  A 
filter).  The  optical  depth  t is  given  by  the  expression 
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Step  No.  3 


Variable  Entering  3 
Flevel  3607.594818 
Standard  Error  of  Y 
Constant 


Variable 

X 1 
X 2 
X 3 

Step  No.  4 


9.500789 

Coefficient 

-1.48664 

-.04167 

.26491 


.001168 

Std  Error 
of  Coeff 

.00210 
.00013 
.00441 


Variable  Entering  4 
Flevel  5.645376 

Standard  Error  of  Y .001157  (oy) 
Constant  9.500496  (a0) 

Std  Error 

Variable  Coefficient  of  Coeff 


X 1 -1.48664  (aj)  .00208 
X 2 -.04058  (b,)  .00048 
X 3 .26491  ( a ) .00437 
X 4 -.00061  (b‘)  .00026 


FIGURE  14.  A Representation  of  the  Computer  Printout 
Showing  Results  for  the  Last  Two  Iterations  in  the 
Determination  of  the  Coefficients  of  the  Polynomial 
PjCx.m)  = aD  + ajX  + bim  + a2X  + b2in2.  The  labels 
in  parentheses  are  not  part  of  the  printout  but  were 
added  to  the  representation  to  assist  in  its  inter- 
pretation. 


2.  Plots  of  the  effective  absorption  coefficient  nij(x,m)  and  its 
polynomial  approximation  Pj (x,m)  against  the  amount  of  ozone  x for  the 
m = 0 case.  Figure  15  is  an  example  of  the  second  kind  of  plot. 
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FIGURE  15.  Effective  Absorption  Coefficient  for  the 
J = 1 Filter  of  the  LD-14  Filter  Assembly  for  the  Case 
id  = 0.  The  calculated  values  of  aj(x,0)  are  given  by 
dots.  The  corresponding  values  for  the  polynomial  are 
represented  by  circles. 


VII.  SUMMARY  AND  COMMENTS 

The  measurement  techniques  and  computer  programs  described  in  this 
report  have  been  used  to  calculate  the  effective  ozone  absorption  coef- 
ficient <j(x,m)  and  the  coefficients  of  its  polynomial  approximation 
Pj(x,m)  for  each  of  the  four  filters  of  the  seventeen  filter  assemblies. 
Several  diffusers  have  also  been  measured  and  their  transmittances 
calculated.  All  but  one  of  the  filter  assemblies  have  been  flown  in 
payloads  carried  aloft  by  Super  Loki  Dart  rockets;  the  polynomials  have 
been  used  by  NASA  to  calculate  ozone  concentrations  from  the  flight 
data. 
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It  is  anticipated  that  the  programs  will  be  used  to  calculate  the 
polynomial  coefficients  for  20  to  25  additional  filter  assemblies  form- 
ing a more  homogeneous  group  than  the  first  seventeen.  Filter  parameters 
such  as  bandpass  position  (midpoint  between  the  half-power  points)  and 
bandwidth,  will  be  tested  to  see  whether  they  may  be  used  to  calculate 
the  polynomial  coefficients  for  a filter  assembly  provided  that  the 
polynomial  coefficients  are  known  for  one  assembly  of  the  group  or 
batch.  Plots  of  the  constant  term  aQ  of  expression  (10)  against  the 
bandpass  position  have  indicated  that  ac  is  related  to  bandpass  position. 
A similar  plot  for  the  a}  coefficient  did  not  lead  to  any  obvious  rela- 
tionship. 

! 
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Appendix  A 

TAPE  PUNCH  CODE  AND  RECORD  FORMATS 


TABLE  A-l.  Tape  Punch  Code  (IBM  Odd  Parity  Code). 


Channel  number/ identity 


8 

7 

6 

5 

4 

o 

3 

2 

1 

Character 

EL 

X 

0 

CH 

8 

o 

4 

2 

1 

0 

X 

o 

1 

o 

X 

2 

o 

X 

3 

X 

o 

X 

X 

4 

o 

X 

5 

X 

o 

X 

X 

6 

X 

o 

X 

X 

7 

o 

X 

X 

X 

8 

X 

o 

9 

X 

X 

o 

X 

# 

X 

o 

X 

X 

- 

X 

o 

CR 

X 

X 

X 

o 

X 

X 

SPACE 

X 

o 

EOA 

X 

o 

x 


= Hole  in  tape 
= Feed  hole 
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TABLE  A-2.  Identification  Number  Record  Format 


Data  recorded 


Identification  number 

(parameter 
switch  settings) 


TABLE  A-3.  Transmittance-Wavelength  Data  Record  Format 


- \ 


Frame 


6 

7 

8 
9 
0 


11 


Data  recorded 


Space 


//  (Ordinate*3  values  above  1000)  or 
- (Ordinatea  values  below  zero)  or 
Hundreds  decade  ordinate  value  (0-9) 


Tens 

Units 


Ordinate3  value 


Space 


Wavelength 

(abscissa) 

value 


CR 


a Transmittance  without  regard  to  the  decimal  point, 
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Appendix  B 

NAMES  OF  ELEMENTS  IN  PROGRAM  FILES 


1.  TAPECARD  FILE 


The  TAPECARD  program  is  filed  under  ROCOZ*TAPECARD  with  the  follow- 
ing element  names: 


FOR  PCH 
FOR  PCHFOR 
REL  PCHFOR 
ABS  PCH 
MAP  SYM 

2.  EDIT  FILE 


This  main  symbolic  program  is  coded  in  FLECS. 
The  FORTRAN  version  of  the  above  program,  PCH. 
The  relative  program  for  PCH. 

The  absolute  PCH  program. 

The  symbolic  input  for  the  MAP  processor. 


The  EDIT  program  is  filed  under  ROCOZ*EDIT  and  contains  the  follow- 
ing elements: 


MAIN  - The  main  program  orders  events,  calls  subroutines,  reads, 
and  writes. 

SETP  - This  subroutine  sets  the  parameters  for  plotting. 

The  system  subroutine  PLT*PLOTS$LIBRARY  must  be  included  in  the  MAP 
to  provide  for  plotting. 

3.  OFFSET  FILE 


The  OFFSET  program  is  filed  under  ROCOZ*OFFSET  with  the  following 
e 1 emen  t name  s : 


MAIN 

EDIT 

AVR1J 

ORDER 


PLOT 

SETP 

MAP 


The  main  symbolic  program  orders  events,  reads  input, 
and  prints  and  punches  results. 

Edits  the  data. 

Performs  a 21-point  running  average. 

Reorders  the  data  so  that  the  first  data  point  consists 
of  the  wavelength  nearest  2500  A,  along  with  the 
corresponding  transmittance  data. 

Plots  the  data. 

Sets  parameters  used  in  plotting. 

Contains  the  names  of  files  and  elements  which  must  be 
collected  to  run  the  program. 
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] 

) 


, 1 


4.  DIFFUSER  FILE 

The  DIFFUSER  program  is  filed  under  ROCOZ*DIFFUSER  and  contains 
elements  having  the  same  names  as  those  in  the  OFFSET  file. 

5.  ALPHAEFF  FILE 

The  ALPHAEFF  program  is  filed  under  ROCOZ*ALPHAEFF  and  contains  the 
following  elements: 


MAIN 


MRL 

*RESTEM 

PLOT IN 

10(A) 

Q(A) 

S (A) 

ALPHA (> ) 
BETA(A) 
MAP (A) 


The  main  program  which  calls  the  subroutines,  •'eads  and 
writes.  It  calls  PLOTIN  for  plotting  the  ALPHAEFF 
function  and  the  polynomial  function  against  x for 
m = mQ,  J = 0,  1,  2,  3.  Additionally,  it  calls  PLOTIN 
to  plot  the  numerator  of  the  integrand  of  the  ALPHAEFF 
function  against  A for  m = mQ , OPDEP  = 0,  1,  2,  3,  4, 
and  J = 1,  2,  3,  where  OPDEP  is  "optical  depth"  defined 
as  aeff(x>m)  multiplied  by  x. 

A subroutine  which  calls  the  RESTEM  subroutine,  furnishes 
the  parameters,  and  prints  the  results. 

A subroutine  which  performs  a multiple-linear-regression 
on  values  furnished  by  ^eff(x>m)  to  fit  a polynomial 
P = aQ  + a i x + b jm  + a2X'  + bpm  . RESTEM  is  called 
by  MLR. 

A subroutine  which  plots.  It  is  called  either  by  MAIN 
or  by  TRY. 

A data  element  giving  extraterrestrial  solar  irradiance 
as  a function  of  wavelength  A,  A = 2500  to  3399  X. 

A data  element  for  photodetector  response  for  900  values 
of  A,  A = 2500  to  3399  X.  This  element  is  replaced 
whenever  the  photodetector  characteristic  is  changed. 

A data  element  providing  transmittance  values  for  the 
diffuser  for  900  values  of  *,  A = 2500  to  3399  A. 

This  element  is  replaced  whenever  the  type  of  diffuser 
is  changed. 

A data  element  providing  absorption  coefficients  for  900 
values  of  A,  A = 2500  to  3399  X. 

A data  element  providing  Rayleigh  extinction  coefficients 
for  900  values  of  A,  A = 2500  to  3399  A. 

A symbolic  element  giving  the  names  of  the  files  which 
need  to  be  collected  for  this  program. 


* RESTEM  is  a subroutine  from  UNIVAC's  STAT-FACK,7  UP-7502,  Rev.  2. 
Sperry  Rand  Corporation.  UNFVAC  Large  Scale  Systems,  STAT-PAOK 
Programmer'  Reference,  UP- 7502,  Rev.  1,  November  1973. 
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Appendix  C 


DISCRETE  VALUES  OF  x,  OZONE  AMOUNTS,  AND  m,  AIR  MASS,  USED  IN 
CALCULATION  OF  a(x,m),  THE  EFFECTIVE  ABSORPTION  COEFFICIENT 


X, 

Ozone  amount 

m 

, Air  mass 

Size  of  No.  of 
increment  values 
DELX  used,  NXS 

Size  of 
increment 
DELM 

No . of 
values 
used,  NMS 

Range 

of 

m 

0.05 

12 

0-0.55 

0.2 

10 

0-1.8 

0.02 

24 

0-0.46 

0.2 

10 

00 

r—H 

1 

o 

0.005 

20 

0-0.095 

0.01 

5 

0-0.04 

0.001 

20 

0-0.019 

0.01 

5 

0-0.04 
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1 New  Mexico  State  University,  Las  Cruces,  NM  (Physical  Sciences 
Laboratory,  B.  G.  Gammill) 

1 Princeton  University,  Forrestal  Campus  Library,  Princeton,  NJ 
1 Rocketdyne,  Canoga  Park,  CA  (Technical  Library) 

1 Rocketdyne,  McGregor,  TX 
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Division) 

Technical  Library 

1 Thiokol  Chemical  Corporation,  Bristol,  PA  (Technical  Library) 

1 Thiokol  Chemical  Corporation,  Elkton,  MD  (Technical  Library) 

1 Thiokol  Chemical  Corporation,  Huntsville  Division,  Huntsville, 

AL  (Technical  Library) 
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(Technical  library) 
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CA  (Technical  Information  Division) 
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